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Abstract 

The third order correlation function of the scalar field advected by a Gaussian random 
velocity, with a spatial scaling exponent 2 — e, and in the presence of a mean gradient, is 
calculated perturbatively in e << 1. This expansion corresponds to the regime close to 
Batchelor's advection by linear diffeomorphisms. The scaling exponent is found to be equal 
to 1 in dimensions 2 and 3, up to corrections smaller than 0(e), implying an anomalous 
scaling of the third order correlation function and the persistence of small scale anisotropy. 



PACS numbers : 

47.27Gs. - Isotropic turbulence; homogeneous turbulence. 
47.27Te. - Convection and heat transfer. 
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The investigation of the statistics of the passive scalar field advected by random flow 
is interesting for the insight it offers into the origin of intermittency and anomalous scaling 
of turbulent fluctuations. The problem studied in this paper is stated simply by: 

d t Q + (u- V)6 = kV 2 G (1) 

with the scalar field O forced by the externally imposed gradient, g. It is convenient to 
subtract out the gradient and study the fluctuating field, 6{r) = 0(r) — gr. It turns 
out that even a Gaussian random, but scale invariant, velocity field results in non-trivial 
anomalous scaling of the passive scalar structure function, < (6(r) — 0(O)) n > for n > 2. 
This has been argued by Kraichnan 1 , via a plausible closure scheme, for a "white noise" 
model where : 

< v a (r, t)v b (r', if) >= 6(t - t')C ab (r - r') (2a) 

with : 

D ab (r) = C ab (0) - C ab (r) = D ((d- 1 + ( v )5 ab - (26) 

(where ( v is the scaling exponent and d the space dimension) which he has introduced some 
30 years ago. 2 The existence of an anomalous scaling has been demonstrated explicitely by 
Gawedski and Kupiainen 3 , and Chertkov et al. 4 for certain limits of this Kraichnan model 
and by Shraiman and Siggia 5 for a generalized phenomenological model where temporal 
correlation of the advecting field is set by eddy turnover. These calculations are based 
on the so called Hopf equations - the stationarity conditions of the equal-time multipoint 
correlators. For the white velocity case these can be derived exactly 6 ' 7 extending the 
original analysis of the 2-point function by Kraichnan 8 . They have the form: 
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(with implicit summation over repeating indices). We restrict ourselves to the inertial 
range of scales, where r is large enough so that the molecular diffusivity can be neglected: 
r»rj = (k/D ) 1/Cv . 

The analysis of Ref.3 is based on the expansion of Eq.(2) in ( v « 1 about the 
diffusion limit ( v = 0, while we consider the complementary limit of ( v = 2 — e, e << 1. 
Reality for the white velocity model, ( v = 4/3 , lies inbetween. The expansion in small e 
is more involved than what was required in Refs 3,4 for two reasons. There are an infinite 
number of degenerate modes for e = which are all mixed by the perturbation, which 
itself is singular 5 ' 9 . That is the perturbation is formally small because of e, but in certain 
restricted regions of configuration space it is the biggest term in the equation. It must 
be treated by the method of matched asymptotic expansions. The exponent we find for 
the third order correlator 10,11 A 3 1 implies that the anisotropy introduced by the mean 
gradient, g on the large scales, decays more slowly as one descends in scale than that 
predicted by K41 theory 12 ' 13 (which for ( v = 2 — e predicts an exponent 1 + e). Since 
the experimental exponent is also approximately one 14 , it will be of interest to compare 
also the full coordinate dependence of the three point correlation function when the latter 
becomes available from experiment or simulations. One way of expressing our result for 
this correlation function, is as an expansion in the degenerate modes of the e = problem. 
Our matching determines all the coefficients explicitly. 

Determination of the anomalous exponents reduces to finding the zero modes of the 
linear operator entering the Hopf equation, 3 ' 4 ' 5 i.e. the left hand side of Eq.(2) which 
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in the present model (and k — > limit) is the generalized Richardson diffusion operator: 
L(dXv) = Y^i^j^abiji — r j)9r.d^,.. The £„ = 2 case is the Batchelor limit 15 which is 
constrained by an overall SL(2) x SO(d) symmetry so that the spectrum of L = L{d, 2), 
also refered to here as the Batchelor-Kraichnan operator, can be completely constructed 
with the help of Lie algebraic methods 5 ' 16 . Here it will serve as a starting point for the 
perturbation theory for £„ = 2 — e: 

L(d,2-e) =L (d) - eL x {d) (3) 

to the leading order in e << 1. (Note, that we are ultimately interested in the physical 
case of Q v = 4/3.) The perturbation expansion around £„ = 2 is a singular problem which 
however can be addressed by the method introduced in Ref.9, as we explain now. 

Let us start with L . It is convenient to introduce the variables p\ = (r*i — f 2 ) / V2 and 
f>2 = (ri + f 2 — 2F3) / v^6 (On occasion we shall refer to i = 1, 2-index labelling the p vectors 
as the "pseudo-space" index to distinguish it from the d-dimensional real space.) Next we 
"factorize": pi = ^ Rii>(x)£,i'Vi> 5 where R represents pseudo-space rotations by x, and fj ly2 
are two orthogonal unit vectors. In d = 3, we also define 7)3 = f]\ A 7)2 = pi A P2/IP1 A p*2 1 
each component of which is invariant under the action of <SX(2) 9 . Another important 
invariant is the area of the r*i, r*2, ^3 triangle: £ = \p\ A p*2 1 = ^1^2- 

The zero modes of Lo(<i) for d = 2, 3 have been constructed in Ref. 5,9; e.g. in d = 3, 
the complete set of eigenf unctions has the form: 

where £ = (£1 + £|)/2£i£2, D l m,m' (v) is the matrix element of the representation of the 
.SO (3) group 17 of order / and PJ' m is the Jacobi function. 18 (Note that quantum number 
m' corresponds to rotations of the r) triad about r/3 in pseudo-space.) To ensure analyticity 
in the £ = \pi A p 2 \ ^0 limit (which corresponds to all 3 points of the correlator being 
on one line), A/2 — max{v, —v — 1) must be a positive integer. This is because as £ — > 0, 
£ ~ £" x -> 00 and PJ' m (£) ~ ^^(".-"-i). 
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The 3-order structure function, or the skewness, which is the physical object of interest 
has odd spatial parity and hence is only non-zero in as much as the mean scalar gradient, 
g, introduces a particular direction. Hence the relevant eigenfunctions are the p-waves, 
1 = 1. The zero mode of L (d) corresponding to the smallest exponent A is obtained for 
\ = v yielding, in the / = 1 sector, A = 1 independent of d. 

We shall need the explicit form of the Lq operator: 

-Lo</tiKx,?7) = %((C 2 - l)d&) + x 4( 3 e ; _^ X il)-v{v + 1)V> (5), 



where v(v + l) = + A) - #rZ(Z + d-2) and I 3 = ^ (r]id 2 - r/ 2 di) . In agreement with 



Eq.(5), the d x and if are diagonalized by exp(zgx) and f]\ ± ir\ 2 the latter corresponding 
to the I = 1, m! = ±1 sector. Requiring the left hand side of Eq.(5) to vanish would make 
it into a Legendre equation 18 with v and hence A entering as an eigenvalue. 
Next we define the perturbation operator in Eq.(3): 

Ll ^ =C + + d ~ 2) + A ' ~ dX ~ d^I L °) (6a) ' 



with 



£=£-ln(| Pl |)x 

(d + i) P i(d^ - U^) - 2 P a lP \{d a 1 d b 1 - l -d a 2 d b 2 \ 



(66). 



In Eq.(6b), the summation extends over all the cyclic permutations of (fi, f 2 , r^), resulting 
in the following symmetry for the p : p\ — > — ^- ± ^p 2 and p 2 — > T^Pi — if- The 
expansion breaks down for — In \p\\ « e . When £ — > oo, i.e., when the three points r*i, 
r2 and r*3 are almost aligned, the operator eC becomes much larger than the Batchelor- 
Kraichnan operator L . This can be seen by expanding the full operator in the limit 
£ — > oo. Defining V'fLi = C^vi^i Vi x)> one finds, for d=3 and the / = 1 sector: 
= 

- ln[l - (1 - r 2 ) cos(2 X )] x [C 2 ^ + £C 1( p + C <p + 0(1/0} (7a). 
+ (X -> X + y) + (X -> X - y) 



The operators £; are : 

C 2¥ > = ^(cos 2 X - l)(4cos 2 X - 1) (-(2^ c - A) 2 + 4 m d m )<P (76), 
£iV?= y cos x sin x(cos 2 x-l) (^(2^-1) (772^ - 7/1^) - 2r] 2 d m ^(p (7c), 



and 



£ V? = " jj( cos2 X - l)(4cos 2 x - 1) {d x 2 + 5)ip 

+ y sin 3 x cos xd x p + ^ (1 + 4 cos 2 x - 8 cos 4 x)Vid m <P 



(Id). 



In Eq.(7d), A has been set equal to 1 - its unperturbed value - since the corrections would 
be higher order in e. The singular nature of the perturbation follows from the fact that the 
C 2 term enters with prefactor £ 2 so that when £ >> (1/e) 1 / 2 , e£ >> L . This situation 
calls for the "boundary layer" type matched asymptotic analysis which we outline below. 

Let us assume A = 1 + e5, define the rescaled 'inner' variable z = e 1 / 2 ^, and introduce 
the function cp(z, x, v) = z ^ (^lfo x)Vi +^2 (-2, x)Vz)- The prefactor is chosen to offset the 
scaling factor £ A / 2 (see Eq.(4)) which vanishes for collinear points. Physics requires that 
<pi is bounded when z — > 00. With this change of variable and functions, the problem can 
be written, provided x ~ £ _1 as : 

(((^ 2 d 2 + 3^) + ^ 2 ^(x)(^^ (8a), 
with 

U( X ) =((cos 2 x - l)(4cos 2 x - 1) Ml - cos(2x)) 

V , (86). 

+(X ^ X + 2tt/3) + (x - X - 2tt/3) J 

The operators L\ and L 2 can be deduced from a systematic expansion of the operator in 
powers of e starting from Eq.(7a-d). 

The boundary conditions at infinity imply that when z — > 00, the solution is a function 
of x only. By direct substitution, one finds that <p\ = and <p 2 = a(x), where a(x) is an 
unknown function, decomposed for convenience as a Fourier series in x : a {x) = ct q e iqx . 
When z — > 0, the problem reduces to the unperturbed Batchelor-Kraichnan operator, up 



to small corrections. In the matching region, defined by z — > but £ — > oo, or equivalently, 
1 << £ << e -1 / 2 , the £ dependence of each Fourier mode in x, q, must match with 
the asymptotic behavior of the eigenmodes of Lq (Eq.(4,5)) which is best found via their 
integral representation given in Ref.5,9,17. One finds that the functions 4>i >q must behave 
as: 4>i^ q = e 2 J g| + 4 ~a + ... and <p2,q = (1 — + ...)sign(q). The crossover equation, 
Eq.(8), can be solved analytically, to the leading order in e, and the imposition of the 
matching conditions determines a(x) via an eigenvalue equation for 5: 

U(x)(d 2 x + l)a( X ) + 65a = (9). 

The analysis in d = 2 can be carried out in a completely similar way. As in the 
3-dimensional case, the behavior for £ — > oo is of the form = a2d(x)£ A ^ 2 + 0(£ A / 2_1 ), 
and the function (i2d is determined by a matching condition. Surprisingly, the equation 
determining a,2d{x)i an d the correction to the scaling exponent, 8 = (A — l)/e is identical 
to Eq.(9). 

Before solving Eq.(9), one needs to determine the appropriate boundary conditions. 
Because the 3-point correlation function must be odd under pi — > —pi, implying that 
a(x + 7r) = — a(x). This, together with the periodicity a(x + 2tt/3) — a(x), resulting from 
the invariance under cyclic permutation of fi, r 2 and f^, implies that a(x + 7r/3) = —a(x)- 
The limit x ^ corresponds to the case where r*i and f 2 come close together : |r 2 — r\ | << 
1^3 — f 5 ! j , |r*2 — ri | . In this limit, the correlation function must be invariant when f± and 
?2 are permuted, implying that a(x) must be even. Since a is even near x = and a is 
antiperiodic with period 7r/3, a(n/6) = 0. 

At small, but finite e, Eq.(9) reduces for x — ► 0, to : 

-X 2 ln(x)a"(x)+5a(x) = (10a). 

Introducing the change of variables : y = — ln(x) and /(y) = a(x), Eq.(lOa) reduces to 
the following (Kummer) equation : 

+ /') + <*/ = o (106). 
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The behavior of the solution when x — * is: 

f(y) ~ y~ s when y — > oo (10c). 

This function diverges (goes to zero) when 6 < (6 > 0). 

Since % — > (for £ _1 = 0) corresponds to p\ — > 0, the perturbation expansion leading 
to Eq.(9) is valid only for y = — lnx << 1/e. Hence, to determine the correct boundary 
condition as y — > oo the solution of Eq.(9) must be matched with the "inner" solution 
describing the correlator with two points near coincidence. The latter is governed by the 
equation derived directly from Eq.(3) by expanding in p\j pi « 1 instead of e and which is 
written conveniently in the polar coordinates |pi/p2| 2 = £~ 2 + X 2 /4 and 9 = arctan(2/£x) 
(restricting here to d = 2 for simplicity). The natural radial variable in this "inner" 
equation turns out to be Y = \pi\ e . The region of matching with Eq.(lOa) corresponds to 
1 — Y « 1 and 6 = 0. Quite generally the solution near Y = 1 behaves as A + B(1 — Y) a 
with a > required to keep the solution from diverging. In the matching region Y 
1 + eln(x/2) so that only the constant term, A, must be kept when computing to the 
leading order in e. Comparing with Eq.(lOc) one concludes that matching the "inner" 
solution is only possible for 6 = 0. 

For 6 = the solution of Eq. (9) is a(x) = sin(-7r/6 — %)/ sin(7r/6) for < x < tt/3 which 
is continued over the full range of x using reflection symmetry and periodicity defined 
above. One observes that a(x) has an apparent |x| singularity near x = (and other 
points related by symmetry) which is regularized only for x < e _1 / e via the crossover 
to the "inner" solution for nearly coincident points pi/p2 < e -1 / e as discussed above. 
Note that although to 0(e) there is no correction to the A = 1 eigenvalue, the computed 
eigenfunction is non-trivial: it is a superposition of many i\)\, q modes since a q ~ 1/q for 
large q. Also note that the calculations in 2 and 3-dimensions are identical and give the 
same result: 6 = 0. 

Thus, the main result of this paper : the scaling exponent of the n = 3 structure 
function behaves as A = 1, up to corrections smaller than e 1 . The exponent of the 3-point 
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correlation function is therefore smaller than the "naive" scaling exponent, equal to 1 + e, 
therefore demonstrating that the behavior of the skewness near the Batchelor limit of the 
Kraichnan, white velocity, model is anomalous. Dispersion in the presence of a mean 
gradient has been shown experimentally 19 ' 20 and numerically 21 ' 22 to give rise to strong 
intermittency effects, resulting in a skewness which remains of order 1, independent of the 
Reynolds number. As it is the case in real flows, it is interesting to notice that even for a 
white noise velocity field, the anisotropy induced at large scales decays more slowly than 
predicted by standard phenomenological arguments. 12 ' 13 ' 23 We conclude by mentioning 
numerical results demonstrating that a large scale anisotropy, such as a large scale shear, 
imposed on a turbulent velocity field, may also result in a large anisotropy at small scale 24 , 
suggesting also the existence of an anomalous exponent for the n = 3 structure function 
of the velocity field. 
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